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DETECTION OF EDGES IN SPECTRAL DATA II. 
NONLINEAR ENHANCEMENT* 

ANNE GELBt AND EITAN TADMOR* 


Abstract. We discuss a general framework for recovering edges in piecewise smooth functions 
with finitely many jump discontinuities, where [f](x) := f(x-\~) — f(x—) 7 ^ 0. Our approach is based 
on two main aspects— localization using appropriate concentration kernels and separation of scales 
by nonlinear enhancement. 

To detect such edges, one employs concentration kernels, K e (-), depending on the small scale e. 
It is shown that odd kernels, properly scaled, and admissible (in the sense of having small W~ 1,0 °- 
moments of order (9(e)) satisfy K e * f{x) = [/](#) + (9(e), thus recovering both the location and 
amplitudes of all edges. As an example we consider general concentration kernels of the form K^(t) = 
52 cr(k/N) sin kt to detect edges from the first 1/e = N spectral modes of piecewise smooth /’s. Here 
we improve in generality and simplicity over our previous study in [A. Gelb and E. Tadmor, Appl. 
Comput. Harmon. Anal., 7 (1999), pp. 101-135]. Both periodic and nonperiodic spectral projections 
are considered. We identify, in particular, a new family of exponential factors, a exp (-), with superior 
localization properties. 

The other aspect of our edge detection involves a nonlinear enhancement procedure which is 
based on separation of scales between the edges, where K e * f{x) ~ [/](#) 7 ^ 0? and the smooth 
regions where K e * / = (9(e) ~ 0. Numerical examples demonstrate that by coupling concentration 
kernels with nonlinear enhancement one arrives at effective edge detectors. 

Key words, piecewise smoothness, concentration kernels, spectral expansions 

AMS subject classifications. 42A10, 42A50, 65T10 


1. Introduction. We discuss a general framework for recovering edges from the 
spectral projections of piecewise smooth functions. Our approach for edge detec¬ 
tion is based on two fundamental aspects—localization to the neighborhood of the 
edges using appropriate concentration kernels and separation of scales by nonlinear 
enhancement. Both the locations and amplitudes of all edges are recovered. 

Let Swf(x) denote the spectral projection of a piecewise smooth /. Given Sjyf, 
one can accurately reconstruct / away from its discontinuous jumps, e.g., 
sect. 2.1], as well as up to the discontinuities (n]]. In either case, an a priori knowledge 
of the location of the edges and their amplitudes is required. This issue was treated in 
recent literature; consult 0],!, 11,11. In 0] , we unified several types of treatments 
as special cases of appropriate concentration kernels, specifically those discussed in Q 
and Q . Here we improve on these results in both generality and simplicity. To this 
end, let [f](x) := /( x+) — f(x-) denote the local jump function and let us consider a 
concentration kernel A' e (-), depending on a small scale e. It is shown that odd kernels, 
properly scaled, an d admissible (in the sense of having small W _1 ’°°-moments of order 
0(e); see eq. ([D])), recover both the locations and the amplitudes of the jumps so 
that 

(1.1) K e *f(x) = [f\(x)+0(e). 

Thus, K e tends to “concentrate” near the singular support of /. 
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Differentiation of e-supported modifiers is one example for local concentration 
kernels outlined in section 2.2.1 . In section 0 we also address the issue of detecting 


edges in global Fourier projections. Given the first 1/e = IV Fourier modes, we seek 
concentration kernels of the form 


Kn( t) 



sin kt. 


It is shown that if the concentration factors <r(£) = £/ii(£) are normalized so that 
/ /x(£)d£ = 1, then Kff(t) is an admissible concentration kernel, Kfix) —> 
[/](&), and the following error estimate holds: 


k=-N ' ' 

The nonperiodic case is studied in section |j. The analogous results for the Chebyshev 
case (consult Corollary |3^) are written as 


< Const ■ 


log IV 

N ' 


7r-\/l — x 2 

N 


Em 



hT' k {x) - [f](x) 


< Const ■ 


log IV 
N 


The special cases of Fourier concentration factors <r, 
were considered earlier in 0, 0, 0 and (l5[ 


(£) ~ sin a/ and er p (£) = ptf v 
Our general framework motivates 


a new set of C^-exponential concentration factors which yield superior localization 
properties away from the detected edges. 

While (1-1) refers to the asymptotic behavior of the concentration kernel as a 
function of the small parameter e j 0, it is essential to recover the exact locations of the 
edges of / for the accurate reconstruction of/. In section || we discuss another essential 
aspect of edge detection, namely, nonlinear enhancement. To this end, one introduces 
a critical threshold, J cr n, for the amplitude of admissible edges and an enhancement 
exponent, p, to amplify the separation of scales in ( 0 ) between the edges, where 
K e * f(x) ~ [/] (x) yf 0, and the smooth regions where K e * /( x) = O(e) ~ 0. Consider 
the enhanced kernel 


KeAfiW 


K e *f if e p / 2 |/C * f{x)\ P > Jcrit, 
0 otherwise. 


Clearly, with p large enough, one ends up with a sharp edge detector where K f> j[f](x) = 
0 at all but (D(e)-neighborhoods of the jump discontinuities. In this sense, the enhance¬ 
ment procedure actually “pinpoints” the location jump discontinuities, allowing an 
accurate reconstruction of /. The particular case p = 2 corresponds to the quadratic 
filter studied in ]l^). ^2|, in the special context of concentration kernels based on 
localized modifiers. 

2. Edge detection by concentration kernels. 

2.1. Concentration kernels. We want to detect the edges in piecewise smooth 
functions. Assume that /(•) has jump discontinuities of the first kind with well-defined 
one-sided limits, f(x±) = lim x ^ x ± f(x), and let [f]{x) := f(x+) — f(x-) denote the 
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local jump function. By piecewise smoothness we mea 


>0 


( 2 . 1 ) 


F x (t) := /(x + ^ /( * 


t) ~ [/](*) 


G BV[ 0, 6) 


V x. 


In practice one encounters functions f{x) with finitely many jump discontinuities, 
and & requires the differential of f(x) on each side of the discontinuity to have 
boun ded variation. For example, if f'(x±) are well defined (for finitely many jumps), 
then H holds. 

We will detect the edges in such piecewise smooth /’s using smooth concentra¬ 
tion kernels , K(t ) = K e (t), depending on a small parameter e. Such kernels are 
characterized by 


(2.2) K e * f(x) — *[f](x) as e -> 0. 

Thus the support of K e * f{x) tends to “concentrate” near the edges of f(x). One 
recovers both the location of the jump discontinuities as well as their amplitudes. 

To guarantee the concentration property of K f . we seek odd kernels, 


(2.3) K e (-t) = —K e (t), 
which are normalized so that 

(2.4) [ K e (t)dt = -l + 0(e) 

Jt> o 

and which satisfy the main admissibility requirement 


(2.5) 


tK e (t)(j){t)dt 


< Const ■ ejl^llBV’- 


Remarks. 

1. For example, if K e (t) concentrates near the origin so that its first moment 
does not exceed 


( 2 . 6 ) 


J\tK e (t)\dt < Const ■ e, 


then it is clearly admissible in the sense that (2.5) holds. We note that our 
admissibility condition also allows for more general oscillatory kernels, K e (t), 
where (2.6) might fail, ye t (p~5| ) is satisfied due to the cancellation effect of 
the oscillations. Consult (2.17) below. 

2. Observe that the admissibility requirement ( |2.5| ) generalizes both properties 
Vs and T± in the definition of admissible kernel |?| Def. 2.1]. 

We state our main result as follows. 

Theorem 2.1. Consider an odd kernel K e (t), JE1), normalized so that ( |2.4| ) 
holds, and satisfying the admissibility requirement (0T Then the kernel K e (t) sat¬ 
isfies the concentration property ( |2.2| ) for all piecewise smooth f’s, and the following 
error estimate holds: 


(2.7) \K e * f(x) - [f](x)\ < Const ■ e. 

1 Here and below we use BV[a , b] to denote the space of functions with bounded variation , endowed 
with the usual seminorm ||^||sv'[a,b] := J a 
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Proof. Using the fact that K f [t) is odd, we have 

K e *f(x)=- I K e (t)(f(x + t) — f(x — t))dt 
Jt>o 


tK e (t ) 


f(x + t) - fix t) - [f](x) 


lt> 0 


dt-[f]{x) [ K e (t)dt 
Jt> o 


Applying (2.4) yields 


( 2 . 8 ) 


I\ e * f{x) - [f]{x) = - [ tK e (t)F x (t)dt + 0{e). 
Jt> o 


By our assumption in (|2.l| ), F x (t) is BV and it is therefore bounded. Consequently, 
in the particular case that the moment bound (213) holds, the first term on the right 
of (2^) is of order O(e), yielding 


\K e * f(x) - [f](x)\ < Const- J\tK e {t)\dt + 0{e) = 0(e). 

In the general case, F x (t) has bounded variation, and the admissibility requirement 
(2.5) implies that the first term on the right of ( |2.8| ) is of order O(e), and we conclude 


I K e * f(x) - [f](x )| < Const\\F x (t)\\ BV ■ e + 0(e) = O(e). □ 

2.2. Examples of concentration kernels. 

2.2.1. Compactly supported kernels. Our first exa mpl e consists of concen¬ 
tration kernels which “concentrate” near the origin, so that (|2.G| ) holds. We consider 
a standard modifier, <j> e (t) := based on an even, compactly supported bump 

function, <f> £ Cq(— 1,1) with </>(0) = 1. We then set 


(2.9) 


I< e (t) = i </>'(-) = 


Clearly, K e is an odd kernel satisfying the required normalization ( p~i| ) 

f K e (t)dt = — f (j)' f - j dt = —<j)( 0) = —1. 

Jt >o 6 Jt >o \ e / 

In addition, its first moment is of order 


/ \tK e (t)\dt = e |s| ■ |<^'(s)|ds = O(e), 

lt> 0 J 0<s<l 


and hence (J2.6^ holds. Theorem 2.1 then implies the following. 


Corollary 2.2. Consider the odd kernel K e (t) = </>'(f), based on even <f> £ 
Cq(—1,1) with </>(0) = 1. Then K e (t) satisfies the concentration property (2.2), and 
the following error estimate holds: 


( 2 . 10 ) 


</>' * f( x ) = [f](x)+0(e). 
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2.2.2. The conjugate Dirichlet kernel. The conjugate Dirichlet kernel, 


K n (t) 


1 

log IV 


D N (t), 


N 

bxit ) := ^ sin kt, 
fc=l 


is an example of an oscillatory concentration kernel. Clearly, A'.y (f ) is an odd kernel. 
Moreover, the normalization ( [Q| ) holds with e ~ |o{ ] N , 


K]y(t)dt 


it =o 


2 

log TV 


E 

odd k's 


l 

k 


— l + 0(e), e 


1 

logiV 


Finally, summing 


N 

D N [t) = sinfct = 
fc=l 


cos | — cos (IV + i)f 

2 sin | 


we find that the first moment of Kn = — .Djv(t)/ log TV does not exceed 



< Const ■ e, 


1 

log AT’ 


so that the requirement ( |2.6| ) is fulfilled. 

Theorem 2.1 then yields the classical result regarding the concentration of conju¬ 
gate partial sums, i sect. 42], ]^3[ sect. II, Thm. 8.13], 


( 2 . 11 ) 


1 

logiV 


D N f(x) 


[/] (x) + O 



We note in passing that in the case of the Dirichlet conjugate kernel, AW(t) does 
not concentrate near the origin, but instead ( |2.G| ) is fulfilled thanks to its uniformly 
small amplitude of order 0(1/ log N). The error, however, is only of logarithmic order; 
consult (T], sect. 2]. 

2.2.3. Oscillatory kernels: General concentration factors. To accelerate 
the unacceptable logarithmically slow rate of the Dirichlet conjugate kernel in ( |2.11| ), 
we consider a general form of odd concentration kernels 


( 2 . 12 ) 





sin kt , 


based on concentration factors cr(/j) wh ich are yet to be determined. Clearly K^(t) 
is odd. Next, for the normalization (|2.4|) we note that 


K%{t)dt = -2 




<t(x) 


>t =0 


k odd 


dx. 


In fact, the above Riemann’s sum amounts to the midpoint quadrature, so that for 
2^ e C 2 [0,1], one has 


't —o 


K^(t)dt = ~2 


r (w) 




k odd 




(2.13) 
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and thus (2A) holds for normalized concentration factors cr(£), 

r 1 


(2.14) 


€ 


-df = 1. 


Consult [[ij for further refinement concerning the assumed regularity of tr(-) (We note 
that er(-) is rescaled here with an additional factor of — n c ompared to [Q.) 

Finally, we address the admissibility requirement (2.5) (and in particular ( [2.6| )). 
To this end, we proceed along the lines of {t|. Assertion 3.3], utilizing the identity 
(abbreviating £ fc = ^) 


N—2 


2 sm(t/2)K^(t) = ^2 (cr(^fc) - 2cr(^ fc+ i) + cr(^ fc+2 )) 


fc=1 


sin (k + 1 )t 
2 sin(t/2) 


- (a-(l) - cr(^jv-i)) 


t 


sin Nt 
2 sin(f/2) 


0(6) - cr(^i)) 


sint 


2 sin(f/2) 


1 


+ <r(£i) COS - - cr(l) COS + - J t 

(2.15) =: h(t) + I 2 (t) + hit) + hit) - hit). 

This leads to the corresponding decomposition of AT^(t) 

er(l) cos (AT + \)t 


K%r(t) — H-N^t) — 


sin; 


Here, R%it) consists of the first four terms on the right-hand side of ( B >, 
/j(t)/2sin(t/2), and it is easily verified that each one of these terms has a 
small first moment satisfying 


(2.16) 


(and consequently, ( |2.5| ) holds); i.e., 
log TV 


/i=0 


\tRx(t)\dt < Const ■ | |u| |c 2 [o,i] ' 


N 


For example, using the standard bound | sin(A;t)/2sin(f/2)| < minjfc, 1 /t}, the contri¬ 
bution corresponding to the first term, fi(f), does not exceed 


it -o 


h{t) 


2sin(f/2) 


dt < 


max a 


N 2 


N—2 

E 

?:=i 


,1/A 

/ 

/t=0 


lt=l/N 


sin(fcf) 


2sin(f/2) 




Similar estimates hold for the remaining contributions of / 2 , h , and h- In particular, 
since ct(£)/£ is bounded, |cr(l/lV)| < Oil/N), and hence / | tXi (i)/ 2 sin(t/2)|dt = 
Oil/N). 

Finally, the admissibility of the fifth term on the right of ( |2.15| ) is due to standard 
cancellation which guarantees that (|2.5|) holds: 


(2.17) 


r (l) 


cos (N + I)t , 
t - \ 2J <j>{t)dt 


' t=0 


< Const • ]</>]] BV 


It is in this context of spectral concentration kernels that admissib ility re quire s 
the more i ntricate property of cancellation of oscillations. Summarizing ( [2.13| ), ( [2.16 ), 
and (2.17), we obtain as a corollary an improved version of the main result in [7], 
Thm. 3.1] regarding spectral edge detection using concentration kernels, K^{t). In 
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particular, since (t) are N-degree trigonometric polynomials, one detects the edges 
of the piecewise smooth function f{x) directly from its spectral projection Sn{ f) '■= 

E n i p ikx 
k=-N JkZ i 


Kjq* f = k n* S N (f) = iir ^2 s 9 n{k)a f k e' 

k=—N ' ' 

Corollary 2.3. Consider the odd concentration kernel (2.12|) 


N 


K N{t) = ~ a ( jy ) Sin H ’ 


k=1 


r M 

s 


G C 2 [0, 1], 


Assume that er(-) is normalized so that (|2.14) holds: 




df = 1. 


Then K^(t) admits the concentration property (2.2), and the following estimate holds: 

log N 


(2.18) 


\ k n * S N (f) - [f](x)\ = Const ■ 


N 


Rem ark. One can relax the regularity on the concentration factor ct(-) 0. Corol¬ 
lary T3 is a generalization of i Thm. 3.1] ;0 in particular, the error estimate (B is 
valid throughout the interval, including at the location of the jump discontinuities. 

Let us introduce few prototypical examples of concentration factors o(-) for the 
detection of edges from spectral data. In this context we note that other detection 
methods of discontinuities in periodic spectral data can be found in the works of 
Eckhoff ||, [§] and of Mhaskar and Prestin, e.g., [^5) and the references therein . We 
note that our results apply to the nonperiodic expansions as discussed in section 3.2.2| 
below. 

1. Trigonometric factors. We consider concentration factors of the form er(£) = 
cr a (£) '■= sin a£/Si(a) with the proper normalization Si(a) := L (sin Tj/rj)dp. 
The edge detector introduced originally by Banerjee and Geer [1| corresponds 
to a n (£); the general case is found in [0, sect. 3.2], 

2. Polynomial factors. As a first example consider cr(£) = £. In this case, 
Kff * f = fS N (f)', and Corollary ^3] recovers Fejer’s result, |2^, sect. Ill, 
Thm. 9.3], with the following error estimate: 


(2.19) 


^M/)'-[/](*) 


< Const 


log N 
N ' 


This is the first member of a whole family of polynomial concentration factors, 
e.g., [0 sect. 3.4], 


( 2 . 20 ) 


a(Z) = aP(0 :=p£* 


which correlate to concentration kernels satisfying ( |2.3| ), (2.4), and (2.5). 
For odd p’s, K^ * f = (—l)^^ 2 ^pN~ p S^\f); for even p’s, K^ * f = 


2 We note the different rescaling here of cr(-) by an additional factor of — 7r, compared with the 
formulation in [M Thm. 3.1]. 
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(—1 ) p / 2 pTV P H * where H(x) = sgn(k)e lkx . These edge detec¬ 

tors were introduced in Ml and were recently analyzed by Kvernadze in ln3|. 
Corollary |S| yields ^ ^ 


N 

^ J2 sgn{k)\k\ p f k e ikx 

k=-N 


[/](*) 


< Const ■ 


log TV 
TV 


The last error estimate is (essentially) first order. It is sharp. It was noted 
in |t]. sect. 3.4], however, that a p, s with higher p’s lead to faster convergence 
rates at selected interior points, bounded away from the singularities of /. 
This leads us to the next example. 

3. Exponential factors. Polynomial concentration factors (of odd degree) corre¬ 
spond to differentiation in physical space; trigonometric factors correspond 
to divided differences in the physical space —co nsult the original derivation 
in Q. Our main result stated in Corollary 245 provides us with the frame¬ 
work of general concentration kernels which are not necessarily limited to a 
realization in the physical space. In particular, we seek concentration factors, 
er(-), which vanish at £ = 0,1 to any prescribed order, 


< 2 - 21 > -jp® -o, 


j = 0,1,2,..., p- 1. 


The higher p is, the more localized the corresponding concentration kernel, 
iV^(-), becomes. Here is why. 

Evaluating K^(t) at the equidistant points tf = 2n£/N, 


K' 


N 



sin 


2irk£ 
TV ’ 


we observe that Kjj(te)/N coincide with the ^-discrete Fourier coefficient of 
er(-); since er(£) and its first p-derivatives vanish at both ends, £ = 0,1, there 
is a rapid decay of its (discrete) Fourier coefficients, ay < Const ■ i~ p , 

\Kn(U)\ < Const ■ ||oi| C p+i[ 0 ,i] 

Thus, for t away from the origin, Kfj{t) is rapidly decaying for large enough 
TV’s. Moreover, we claim an increasing number of moments of A^(-) vanish. 
To this end we consider the odd moments of TVy(-) (its even moments vanish, 
of course). With £*, = k/N we find 


JV 


( 2 . 22 ) 


Ex [-J2 


k= 1 


(T ( — | sin kt I dt 


2 N 1 /»7T AT 

= W M I rja ^ sin ^r)dT 


N ° ti N Jt=0 


+2 r 1 dt r vN 

= 1 =o a{0 ' W J C0S ^ T)dTd ^ j ° dd 


T =0 
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I ntegr ate by parts; resp ectively, sum by parts the summation on the right of 
([2. 22). Thanks to (2.21) the boundary terms vanish and we have 


t j K^(t)dt ~±j f ^-a(0 • Sm(7 ^V 


1 

w 


/«=o 

f 1 d p ~ j 


S 


(If®) 


Ji=o '’C ' 

j < p, p odd. 


siiifriVg) /J_ 

( t tN)p-3 \NP 


As an example, we consider the exponential concentration factors 


(2.23) a exp (f) = Const ■ £e, Const = / exp 


-1 


ag(v - 1) 


dp 


normalized so that f^ =0 cr exp (f) / fdf = 1. Here, the Cq° [—1,1] concentration 


factor a exp {f) vanishes exponentially at both ends; £ = 0,1 so that (2.21) 
holds for all p’s. Figures pT3 and 2.4 confirm the improved localization of 


these exponential concentration factors. 

4. Band pass filter. Bauer || has considered a family of what he termed as 
“band pass filter,” supported in the range of middle frequencies, say, 

suppp C [1/4, 3/4]. We note in passi ng tha t these are special cases of p-order 
admissible concentration factors (see ( ED ). although the normalization used 
in eq. (1.35)], f r/(x)/x sin(wx)dx = 1, prevented the recovery of the 
amplitude of the jumps. 

To demonstrate the detection of edges by the concentration factors outlined above, 
consider the following two examples of discontinuous /’s (defined on [— 7r, 7r]): 


fa{x) := -sgnx -cos ^(2 + sgnx)j , f b (x) := 


cos(x — ^sgn(\x\ — ^)), x < 0, 
cos(|a; + xsgn(\x\ — f)), x > 0. 


In both cases, f a (x) and f b (x) are recovered from their Fourier coefficients using the 
Fourier partial sums •Sjv [/](*), and we wish to recover their jump discontinuities 


—2, x = 0, 


[/«](*) = 


0 


else, 


mx) = 


±\/2, x = ±f, 


0 


else. 


Figures 2.1 and 2.2 demonstrate the use of trigonometric and polynomial concen¬ 
tration factors for the detection of edges from Fourier spectral data. 

As noted in sect. 3.4], polynomial factors of higher degree yield improved 
results away from the jump discontinuities. Indeed, the corresponding concentration 
kernels, I\fj (•), have additional vanishing moments. In the limit, one arrives at 
exponential factors, (•); Figures ^ and demonstrate the edge detection of 

these factors in Fourier expansions Sisrfa(x) and 5jv/&(a;)- The improved localization 
is evident, due to the faster convergence rate in the smooth parts of /’s. In particular, 
the superiority of the exponential factors is illustrated in the figures on the left, when 
compared with the first-order accurate polynomial concentration factor, cr p=1 (£) = 
At the same time, Gibbs oscillations can be noticed in the vicinity of the jump 
discontinuities. 
















1398 


ANNE GELB AND EITAN TADMOR 


[f](x) WM 




Fig. 2.1. Trigonometric concentration factor <j(£) = f or fa{x), where the exact 

jump value is [/](0) = — 2 and (right) fb(x), where the exact jump values are [/](d=^) = d=v2. 


MM M(x) 




Fig. 2.2. Jump value obtained by the polynomial concentration factor a p 1 (^) = £ for (left) 
f a (x) and (right) fb(x). 



-2 


'. 

-3-2-10 1 2 3 



Fig. 2.3. Edge detection using the exponential concentration factor a (£) = 3exp(^^—jy) (left) 
a exp vs. a p=1 for S 4 of a (x) and (right) cr exp for Sjv/a(aO with N = 20,40,80 modes. 
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Fig. 2.4. Edge detection using the exponential concentration factor cr(£) = 3exp( ) | , ) (left) 
a ex P vs. cr p=1 for S.]ofh (x) and (right) a exp for S.\ftf.x) with N = 20,40,80 modes. 


3. Edge detection in nonperiodic projections. Consider a piecewise smooth 
/(•). To simplify our presentation, we assume / experiences a single jump discontinu¬ 
ity at x = c. The localization property of the appropriate concentration kernel in the 
presence of a single jump applies to the case with finitely many jump discontinuities. 
We begin with an alternative derivation of our results for the periodic case. 

3.1. Revisiting the periodic case. If a 27r-periodic /(■) experiences a single 
jump, [/](c), then it dictates the Fourier coefficients decay H, [p3[ , 


fk = [/](c) 


„ — ikc 


27 rik 




To extract information about the location of the jump from the phase of the leading 
term, we examine the special concentration kernel, with <r(£) = £, where K ^ * 

s N (/) = §s N (fy, 


p ik(x—c ) 


2irik 




(3.1) = .r £ §{[/](e) 

k=-N ^ 

(3-2) = [/](c)E E { 1 + ° (1) } eiHx ~ C) = W (*) + O fl ° gN 


\ N 

Here we used the concentration property of the Dirichlet kernel localized at x = c: 

1 


— y 

9 IV ^ 


N 


o 


ik(x-c) 


2 N 


k=—N 


iVjx — c| 

1, 


x^c, 


The same property applies to the class of concentration factors, er(£) := £/u(£), such 
that (2.14) holds, f* (i(£)d£ = 1, 


1 

2TV 


N 


E A* 

k—-N 


O 


ik(x-c) _ 


1 


N\x — c\ 

1, 


, x^c, 
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It then follows that the corresponding K^ 
kernel, so that K jy * /( x) —» [f]{x)- 


in (2.12) is an admissible concentration 


3.2. Nonperiodic expansions. 

3.2.1. General Jacobi expansions. We begin with the Jacobi expansion of a 
piecewise smooth /(•), 


(3.3) 


N r l 

S*(/) = £M(*>. fk= f{x)w{x)P k {x)dx. 

fc=0 - ' _1 


Here Pkix) are the Jacobi polynomials—the eigenfunctions of the singular Sturm- 
Liouville problem 


(3.4) 


((1 - x 2 )ui{x)P' k {x))' = -A ku(x)Pk(x), -1 < x < 1 


with corresponding eigenvalues A*, = A^ = k(k + 2a + 1). Different families of 
Jacobi polynomials are associated with different weight functions u>(x) = u> a (x) := 
(1 — x 2 ) a . To simplify the computations, we assume that the Pk s are normalized so 
that ||Pfc(a;)|| w = 1. 

As in the periodic case, integration by parts (against (3T)) shows that a single 
jump discontinuity, [/] (c), dictates the decay of the Jacobi coefficients, 


(3.5) 

A = ^ y i /(*)((! - x 2 )uj(x)P^(x))'dx = [/](c)-^-(l - c 2 Mcmc) + O . 

To extract information about the location of the jump, we consider the conjugate sum 
of the form 


(3.6) 


7T\/l ^ X 2 

N 





fkPk(x) 


= im 


7T\/l — X 2 (l — C 2 )w(c) 

N 


Em 




X P k(c)Pk{x ), 


corresponding to concentration factors er(£) = £/z(£). We shall focus our attention on 
the particular case p(£) = 1, 


(3.7) 


r\/l — ; 

iv - 


-S N (f)'{x) = 


r\/l — . 


N 


N 


■E* p fc( a 


k—l 


= im 


7rVl — x 2 {l — c 2 )uj(c) 


N 


N 


k= 1 




This is the nonperiodic analogue of the Fourier concentration kernel K ^ * / with the 
additional prefactor weight of y/1 — x 2 . 

We want to quantify the localization property of the last summation. To this 
end we note that if {Pf. a \x)} are the Jacobi polynomials w.r.t. the weight function 
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u> a (x), then {P k (x)} are the Jacobi polynomials w.r.t. the modified weight function 
u>p(x) = (1 — x 2 )u a (x) with (3 := a + 1: indeed, their wg-orthogonality follows from 
integration by parts of (3d) against ■ Thus, 


P, 


(a)/ 


(x) = C k<p Pj?\, 


(3 = a + 1. 


The coefficients C/c ,/3 are determined by normalization, where by using (3.4) once more 
we find 


1 = IIP 


(«)||2 


j ^((1 -x 2 )u a {x)pj. ah ypjf‘ ) {x)dx = 


|| p (/3) || 2 
\ " r k— 1 HoJ/3 > 


and hence we set C k ,/3 = so that {Pj_ i} is the orthonormal family w.r.t. up 

weight. Inserted into the leading term of ( p.7| ), we end up with a Jacobi kernel 
associated with weight function u>p(x) = (1 — x 2 ) 13 , 


r\/l — . 

N~ 


-S N (f)'(x) ~ [/](c 


= [/](c 


7T\/l — X 2 0Jf} (c) 

iv 

7rVl — x 2 Ufi (c) 

N 


N 




fc=l 

AT 


?(/^) p(/^) 




/c=l 


We rewrite this as 
7r 


(3.8) 


a/T -- ax 


IV 


^(/)'(x) ~ [/](c 


Tt\Jl — X 2 LOp(c) 


N 


x I\n{c , z). 


By virtue of the Christoffel-Darboux formula, e.g., [|T9], Thm. 3.2.2], the kernel 
Kn(c,x) is given by 


k N -1 Pjp ( x )Pn-i(c) - PX” (c)PriW fcN-1 1 


,(/3) 


(/3)/ 


,(/3) 


(3.9) Kn(c,x) = 


kN 


x — c 


/cat 


and it remains to quantify the concentration property of Kn{c,x). To this end we 
use the asymptotic behavior of P^ which is stated as| 


(3.10) P^\x) ~ min 


KU}p + i/ 2 {x) ’ 


Const ■ N 0+1 > 2 = 


nup +1/2 (x N ) 


where Xn '■= sgn(x ) • min{|x|,l — Const/N} denotes the separation bet ween the 
interior and boundary regions. Using this to upper bound Kn{c,x) in (T9), we find 


(3.11) 


\K n (c,x)\ < 


7r V W /3+!/2 ( C N )w/3+1/ 2 {xn ) \x-C\ 


X ^ c. 


3 The first term on the right of ( J3.10| ) follows from the classical asymptotic formula, e.g., ||l£| 
Thm. 12.1.4], which tells us the behavior of the Lj,-normalized P^\x) at the interior x = cos 6 


yJup +1 / 2 (x)P < ^\x) ~ yj 2/tt cos{N6 + 7 ( 6 )}, x = cos 9. 

The second term on the right of (3.10) reflects the fact that as x approaches the ± 1-boundaries, the 
L^j -normalized P^\x) approaches its maximal value, e.g., Jrj, eqs. (4.7.3), (4.7.15)] 


pg>W = ( N+ %- 1 


r(N + 1)(IV + A) 
7r2 1 - 2A r(Ar + 2A) 


T(A) 


T(A + 2A)(Ar + A) 


JV! 


~N X , \ = 13+ 1/2. 
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The upper bound on the right is in fact the leading term in the asymptotics of Kn(c, x ) 
for large N' s as long as — 1 < x ^ c < 1, e.g., [19, sect. 13.4]. Similarly, the behavior 
at x = c yields 


(3.12) 


K n (c, c) 


N + Const 


-1 < c < 1. 


7rw /3+l/2 (c) 

The desired concentration property now follows, similar to the localization of the 
periodic Dirichlet kernel D^{x — c)/N in (3.3). We restrict our attention to interior 
jumps, —1 < c < 1, so that for N large enough, = c, and ( |3.11 ), ( 3.1 2| ) , and (3.8) 
yield 


(3.13) 

7T\/l — X 2 U 10 (c) 

N 


O 


x Kn{c, x) 


V / <^q+l/2(c) ^ 1 

\J w a+l /2 (Xn) N\x — c| 

1, 


jya+1/2 

N 


, i/c, 


X = c. 


We summarize by stating the following. 


Corollary 3.1. Let S^r)/) denote the truncated Jacobi expansion (T3) of a 
piecewise smooth f, associated with a weight function u> a = (1 — x 2 ) a , — 1 < a < 0. 
Then 7rVT — x 2 SN(f)'(x)/N admits the concentration property 


(3.14) 


7T-\/l — X 2 


N 


SwifYix) - [/]( x) 


< Const ■ 


logN 


N(l — a: 2 ) a / 2 + 1 /4 ’ 


— 1 + Const ■ —p < x < 1 — Const ■ ——r. 

N- TV z 


It is instructive to examine the above discussion for the special case of Chebyshev 
expansion corresponding to a = — 


N or 1 

S N (f)(x)=Y / 'hT k (x), f k = - 


f(x)T k (x ) 
=-i Vl - x 2 


dx. 


(Observe that except for Chebyshev expansion, the concentration bound ( 3.13| ) dete¬ 
riorates as we approach the boundaries, depending on whether |rc| ~ 1 for a > — ^ or 
|c| ~ 1 for a < — i.) The conjugate sum corresponding to (T8) reads 


T\Jl — x 2 

N 


_ N o ,_ ,_ 

SnU)\x) ~ [f](c)-Y,'^VT^T'(c)Vl^T'(x). 
V k =1 


In this case, we can sum the corresponding Chebyshev kernel: setting x = cos 9 and 
— 1 < c = cos rj < 1 we find 


(3.15) 


T\/l — : 

N~ 


N 


SN(f)'(x) ~ [f](c)—^2' am(krj)sm(kd) 


k—0 


= [/](c)^p. n{0-v) -D N {0 + rj) 


O(jf), X^c, 
[/](c), X = c. 
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3.2.2. Chebyshev expansion. Our discussion above on edge detection in the 
nonperiodic expansions is based on expansion of the Jacobi coefficients to their leading 
order in (3.5). More precise information is obtained using the general framework 


introduced in the main theorem, Theorem 2.1 


Corollary 3.2. Let /(•) be a piecewise smooth function with Chebyshev expan¬ 
sion Sn/(x) ~ Y' f(k)Tk(x). Consider the concentration factors, a(f) := £/i(£), with 
fi(-) normalized so that 


/%m = i, M£)eC 2 [ o,i]. 

JO 

Then Kfj(t) * f(cosd) admits the concentration property ( fc.2| ), and the following esti¬ 
mate holds: 


(3.16) 


7T\/l — X 2 

N 


£^ 



hk)T'(x) - mo t) 


= Const ■ 


log N 
N 


Proof. With a piecewise smooth /( x) defined over the interval [—1,1] we utilize 
the usual Chebyshev transformation x = cos 9, 0 < 6 < tt. We consider the even 

extension /(cos0), — 7r < 9 < tt. Using Theorem |2.l| along the lines of Corollary |2.3| , 

we find that the odd concentration kernel, Kfj{f), recovers the jumps of f(cos9), i.e., 


A, U f U \ 

- £ ( N ) sin ( fc ^ * /( cos6) ) “ [/K cos6 0 

fc—1 1 ' ' 


< Const ■ 


logiV 

N 


With T^{x) = — fcsinfc0/-\/T — x 2 , a straightforward computation shows the sum on 
the left equals 


-N / i \ N i i \ 

“ £ ( TV ) sin ( fci ) * f( cos °) = “ 7r £ 'jjV ( Jj ) /(*) sin(fc6») 

k =1 ^ ' k=l 1 ^ ' 


r\/l — ; 

N~ 


N 


■ £ M ( ^ ) f(k) T k(x) 

k=1 V 


and the result follows. □ 

We turn to numerical examples. The following tables summarize our results for 
the edge detection in Legendre expansion, corresponding to a = 0, and in Chebyshev 
expansion, corresponding to a = —1/2. Scaled to the unit interval [—1,1], we consider 
fa(~) and fb( — )- The results confirm the linear convergence rate stated in Corollary 


3.1 


both a way from the jumps—consult Tables 3.1 and |3.2| -as well as at the jump 
itself (Table |ib3|) . 

We note that the critical threshold must be very high for N = 20,40 to eliminate 
the artificial jumps. This indicates that 40 nodes are not enough to resolve the jumps 
of fb{x) in either the Chebyshev or Legendre case. 

It is clear from Tables |3.l| and |3.2| that convergence is nonuniform at the bound- 


aries. We have observed in our numerical experiments that the edge detector, 


N 


S N \f\' 


experiences larger oscillations near the boundaries which do affect the linear conver¬ 
gence rate there. In this context we note the dependence of the error bounds on the 
smoothness of f(V 1 — x 2 ). 


The first-order convergence is reconfirmed, in Table 3.4 below, when measuring 
the L 1 -error away from the jumps discontinuities (and up to the boundaries). 
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Table 3.1 

Pointwise error estimate tt y 7 I - H? /NS N (f a y(x/n) - [fa](x/ 7r) | away from the jump discon¬ 
tinuity at x = 0. 


N 

Legendre expansion 

Chebyshev expansion 


x = —.998 

x = - .9805 

x = .5 

x = .75 

x = —.998 

x = -.9805 

x = .5 

x = .75 

40 

.192 

.112 

.135 

.123 

3.6E-2 

7.1E-2 

.121 

.147 

80 

7.7E-2 

3.8E-2 

.16 

1.3E-2 

2.1E-2 

5.9E-3 

.145 

2.4E-2 

160 

7.6E-3 

4.0E-3 

3.1E-2 

5.6E-3 

5.3E-3 

8.7E-3 

3.3E-2 

1.6E-2 


Table 3.2 

Pointwise error estimate |7r\/l — i?/Ns N {f b y(x / 7r ) — [/b]( aj / 7r ) I awa V from the jump disconti¬ 
nuities at x = 


N 

Legendre expansion 

Chebyshev expansion 


x = —.998 

x = -.9805 

x = .5 

x = .75 

x = —.998 

x = -.9805 

x = .5 

x = .75 

40 

8.0E-2 

7.4E-2 

5.4E-3 

.45 

1.3E-02 

1.5E-02 

9 . 8 E -02 

.56 

80 

.15 

1.55E-2 

7.0E-4 

.18 

1.1E-02 

1.3E-02 

4.5E-02 

.26 

160 

1.6E-3 

7.2E-3 

3.9E-4 

.12 

1.8E-03 

4.6E-03 

2.3E-02 

.12 


Table 3.3 

Pointwise error estimate |7r\/l — x‘ 2 /NSN(f)' (c) — [/](c)| at the point(s) of discontinuity, x = c. 



Legendre 

Chebyshev 

N 

/«(0) 

M-w/2 ) 

M*/2) 

/«(0) 

M-n/2) 

AW 2 ) 

40 

1.9E-02 

8.5E-02 

.187 

IE-2 

5.3E-2 

6.2 E-2 

80 

1.4E-03 

3.9E-02 

.11 

IE-2 

4.7E-2 

2.5E-2 

160 

1.3E-02 

1.7E-02 

5.5E-02 

8.7E-3 

1.1E-2 

2.2E-2 


Table 3.4 

Z/ 1 [— 1,1] — {c} error estimate \ny/l — x 2 /NSN{f)' (x) — [/](»:) | away from discontinuities. 


N 

Legendre 

Chebyshev 


f = fa 

f = fb 

f — fa 

f = fb 

40 

.16 

.25 

.156 

.27 

80 

8.3E-2 

.13 

8.4E-2 

.12 

160 

4.4E-2 

6.8E-2 

4.6E-2 

6.9E-2 


4. Nonlinear enhancement. The detection of edges in Theorem 2.1 is based 
on separation of scales. Thus, consider for example a piecewise smooth / with finitely 
many jump discontinuities at x = ci, C 2 ,..., c n . If K t is an admissible concentration 
kernel, then \K e * /(a;)| -C 1 for x away from these jumps, whereas at x = Cj, K e * 
/(®)*=c ~ 0(1). 


(4.1) 


Ke * f{x)= 


0(e), 

[/] fe), 


X ^ Ci,C 2 , ... ■ 
X = Cj. 


The last statement refers to the asymptotic behavior of the concentration kernel 
as a function of the small parameter e J, 0. In this section we outline a new, nonlin¬ 
ear enhancement procedure, which is easily implemented to “pinpoint” finitely many 
edges in piecewise smooth f’s. 


To this end we enhance the separated scales in (4T) by considering 


(4.2) 


- p HK e *f( X )r = 


~ Const • £ 2 , 
([/]fe')) p -e-| 


X 7^ Ci,C2, ..., 
X Cl 7 C2 5 . . . . 
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By increasing the exponent p > 1, we enhance the separation between the vanishing 
scale at the points of smoothness (of order 0(e ?)) and the growing scale at the jumps 
(of order 

Next, one must introduce a critical threshold which will eliminate all the unac¬ 
ceptable jumps. Only those edges with amplitudes larger than the critical threshold, 
[f](x) > will be detected. Thus J = J cr u is a measure which defines the 

small scale in our computation of edge detection. We note that J = J cr it is data 
dependent and is typically related to the variation of the smooth part of /. 

Given this critical threshold, we form our enhanced concentration kernel K e< j[f] 


(4.3) 


AV [/](*) 


K e *f(x) if e a \K e * f{x)\ v > Jcrit, 
0 otherwise. 


Clearly, with p large enough, one ends up with a sharp edge detector where 

Kt,j[f](x) = 0 at all but 0(e) neighborhoods of the jumps x = ci, C 2 ,_In practical 

applications, a moderate enhancement exponent, p < 5, will suffice. We consider two 
examples. 

1. The quadratic filter. Consider the peaked concentration kernel ( p.9| ) K e (t) = 
fi' e (t). Then, with p = 2, one finds the so-called quadratic filter jl2) , p^ |, 
where 


(4.4) 


{K e * f(x)) 2 = (cj>' e * f(x)) 2 —> \f] 2 (x). 


2. Enhanced spectral concentration kernels. We apply the procedure of nonlin¬ 
ear enhancement in conjunction with spectral concentration kernels K e = 


K^(t) := — XifcLi cr(-^)sinfcf by considering the corresponding enhanced 
spectral concentration kernel 


h' a — 


Kn * A/v fix) 

0 


if e 2 \K$f * S N (x)\ p > Jcrit, 
otherwise. 


The enhanced spectral concentration kernel depends on four ingredients which 
are at our disposal: 

• The number of modes, IV; 

• The enhancement exponent, p; 

• The critical threshold, J; 

• The concentration factor, er(£). 


Figures 4.1 and 4.2 demonstrate the enhancement procedure to the spectral de¬ 


tection of edges depicted earlier in the corresponding Figures and 2.1. 

We conclude with nonperiodic examples. In Figure [h^ we show the detection of 
a single edge in f a (x/n) from its Legendre expansion, J2(fa)kPk{x). The detection 
in Chebyshev expansion is shown in Figure L4 for fb(x/n). In both cases we used an 


enhancement factor p = 2 and a critical threshold J = 5. 

5. Concluding remarks. Accurate reconstruction of piecewise smooth func¬ 
tions from their spectral projections is plausible only when the location (and ampli¬ 
tude) of the underlying jump discontinuities are known; consult §, 0, 0, & 
fl6| , J2l| , ||ll[, and references therein. 

Theorem 2.1, and its Corollaries 2.3 and 3.1, provide the general framework for 


the detection of edges from spectral data, in both periodic and nonperiodic cases. 
The detection is based on admissible concentration kernels which include as particu¬ 
lar cases classical examples of Fejer as well as additional examples in recent literature 
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MM 


0.3 

-0.2 

-0.7 

-1.2 

-1.7 

-2.2 



-0.1 0.9 1.9 2.9 


Fig. 4.1. Jump value obtained by applying the polynomial concentration factor cr(£) = £ with 
N = 40 with enhancement exponent p = 2 for (a) f a {%), where the exact jump value is [/](0) = —2 
and (b) where the exact jump values are [/](±^) = d=\/2. 


m<x) 

0.3 L . J=5 



"-3.1 -2.1 -1.1 -0.1 0.9 1.9 2.9 


-3.1 -2.1 -1.1 -0.1 0.9 1.9 2.9 


Fig. 4.2. Jump value obtained by applying the trigonometric concentration factor cri(£) = 
with N = 40 modes and enhancement exponent p = 2 for (left) f a (x), where the exact jump value 
is [/](0) = —2 and (right) fb(x), where the exact jump values are [/](d=^-) = d=\/2- 




Fig. 4.3. Detection of edges in Legendre expansion of f a ( x/ty) with exact jump value is [f a ]( 0) = 
—2 (left) before and (right) after enhancement with p = 2 and J = 5. 
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Fig. 4.4. Detection of edges in Chebyshev expansion of fi(x/ir) with exact jump value is 
[/i,](±^) = ±\/2 (left) before and (right) after enhancement with p = 2 and J = 5. 


|, In particular, we introduce here a new family of exponential concentra¬ 
tion kernels, (2.23), with a superior convergence rate away from the edges. A linear 
convergence rate is observed near the detected edges. We also introduce a nonlinear 
enhancement (4.3) procedure which enables one to “pinpoint” edges with amplitude 
larger than a critical threshold. 

Recently, the edge detection and enhancement method was applied to nonlinear 
conservation laws ||(J , 0 as a postprocessing tool to improve the overall convergence 
rate of the spectral viscosity solution. Since the edge detection occurs only at the 
postprocessing stage, very little cost is added to the procedure yet the results are 
dramatically improved. Future applications, in both one- and several-space dimen¬ 
sions, will also include image processing, where edge detection is needed to denoise the 
contamination by the C ) (l)-Gibbs oscillations in the neighborhoods of the undetected 
edges. 
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